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Abstract 



The equations of in-medium gluodynamics are proposed. Their classical 
lowest order solution is explicitly shown for a color charge moving with con- 
| stant speed. For nuclear permittivity larger than 1 it describes emission of 

Cherenkov gluons resembling results of classical electrodynamics. The val- 
ues of the real and imaginary parts of the nuclear permittivity are obtained 
from the fits to experimental data on the double-humped structure around 
the away-side jet obtained at RHIC. The dispersion of the nuclear permittiv- 
ity is predicted by comparing the RHIC, SPS and cosmic ray data. This is 
■ important for LHC experiments. Cherenkov gluons may be responsible for 

the asymmetry of dilepton mass spectra near p-meson observed in the SPS 
CN ■ experiment with excess in the low-mass wing of the resonance. This feature 

is predicted to be common for all resonances. The "color rainbow" quantum 
t^j- \ effect might appear according to higher order terms of in-medium QCD if the 

0^ ■ nuclear permittivity depends on color. 

(N 

^ ■ 1 Introduction 

o 



The collective effects observed in ultrarelativistic heavy-ion collisions at SPS and 
RHIC [H El El IU El E] have supported the conjecture of quark-gluon plasma (QGP) 
formed in these processes. The properties and evolution of this medium are widely 
debated. At the simplest level it is assumed to consist of a set of current quarks 
and gluons. It happens however that their interaction is quite strong so that the 
notion of the strongly interacting quark-gluon plasma (sQGP) has been introduced. 
Moreover, this substance reminds an ideal liquid rather than a gas. Whether per- 
turbative quantum chromodynamics (pQCD) is applicable to the description of the 
excitation modes of this matter is doubtful. Correspondingly, the popular theoret- 
ical approaches use either classical solutions of in- vacuum QCD equations at the 
initial stage or hydrodynamics at the final stage of its evolution. However, it is sur- 
prising that no attempts to write down the equations of in-medium QCD, similar to 
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the very successful approach in electrodynamics, were done before the recent paper 
[7] appeared. 

The collective excitation modes of the medium may play a crucial role. One 
of the ways to gain more knowledge about the excitation modes is to consider the 
propagation of relativistic partons through this matter. Phenomenologically their 
impact would be described by the nuclear permittivity of the matter corresponding 
to its response to passing partons. Namely this approach is most successful for 
electro dynamical processes in matter. Therefore it is reasonable to modify QCD 
equations by taking into account collective properties of the quark-gluon medium. 
For the sake of simplicity we consider here the gluodynamics only. The generalization 
to quarks is straightforward. 

The classical lowest order solution of these equations coincides with Abelian 
electro dynamical results up to a trivial color factor. One of the most spectacular of 
them is Cherenkov radiation and its properties. Now, Cherenkov gluons take place 
of Cherenkov photons [SI E]. Their emission in high energy hadronic collisions is 
described by the same formulae but with nuclear permittivity in place of the usual 
one. Actually, one considers them as quasiparticles, i.e. quanta of the medium 
excitations with properties determined by the permittivity. 

These formulae are used for fits of experimental data on the double-humped 
structure around the away-side jet obtained at RHIC. The values of the real and 
imaginary parts of the nuclear permittivity are determined. Comparing the RHIC, 
SPS and cosmic ray data one can guess that the nuclear permittivity depends on 
energy. This leads to predictions for future LHC experiments. 

Beside the high energy region, the effects due to Cherenkov gluons might be no- 
ticed in asymmetry of the shapes of resonances passing through the nuclear medium. 
This prediction is compared with experimental data of SPS. 

Another possible effect of "color rainbow" might appear if the nuclear permit- 
tivity differs for partons of different colors. The higher-order non-linear terms of 
in-medium QCD equations are in charge of it. 

2 Equations of in-medium gluodynamics 

At the beginning let us remind the classical in-vacuum Yang-Mills equations 



(1) 



= d"A v - d u A> x - ig[A", A u ], 



(2) 
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where = A^T a ; A a (A® = $ a ,A a ) are the gauge field (scalar and vector) po- 
tentials, the color matrices T a satisfy the relation [T a , T b ] = if a b c T c , = — 
ig\A^ •], J u (p, j) is a classical source current, h—c—1 and the metric tensor is 

0""=diag(+,---). 

In the covariant gauge d^A^ = they are written as 

HA" = J^ + tg[A u ,d u A^ + F u % (3) 

where □ is the d'Alembertian operator. It was shown [10] (and is confirmed in what 
follows) that in this gauge the classical gluon field is given by the solution of the 
corresponding Abelian problem. 

The chromoelectric and chromomagnetic fields are 

= (4) 
B n = -eWfj; (5) 

or as functions of gauge potentials in vector notations 

dA 

E a = -grad$ a - + gf abc A b $ c , (6) 

B a = curlA a - ^/ a6c [A 6 A c ]. (7) 
The equations of motion (0Q) in vector form are written as 

divE a - gf a6c A 6 E c = p a , (8) 
<9E 

curlB a - - gf abc (<S> b E c + [A 6 B C ]) = j„. (9) 

The Abelian equations of in- vacuum electrodynamics are obtained from Eq. ([3]) if 
the second term in its right-hand side is put equal to zero and color indices omitted. 
The medium is accounted if E is replaced by D = eE in F^ v , i.e. in Eq. (j4j@. 
Therefore the Eqs. (JHJ), (EJ) in vector form are most suitable for their generalization 
to in-medium case. The equations of in-medium electrodynamics differ from in- 
vacuum ones by dielectric permittivity e ^ 1 entering there as 

A. 

AA-e— = -j, (10) 



2 e denotes the dielectric permittivity of the medium. It is well known [11] that magnetic 
properties of a substance are reproduced with the proper account of temporal and spatial dispersion 
of e. 
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<^-^) = ~p- ( n ) 

The permittivity describes the matter response to the induced fields which is as- 
sumed to be linear and constant in Eqs. ([TO]) . (jTTj) . It is determined by the distri- 
bution of internal current sources in the medium. Then external currents are only 
left in the right-hand sides of these equations. 
Now, the Lorentz gauge condition is 

divA + e— = 0. (12) 
at 

The Lorentz invariance is lost if e ^ 1 in front of the second terms in the left- 
hand sides. Then one has to deal within the coordinate system where a substance 
is at rest. The values of e are determined just there. To cancel these requirements 
one must use Minkowski relations between D, E, B, H valid for a moving medium 
[12J. It leads to more complicated formulae, and we do not use them in this paper. 

The most important property of solutions of these equations is that while the 
in- vacuum (e = 1) equations do not admit any radiation processes, it happens for 
e / 1 that there are solutions of these equations with non-zero Poynting vector. 

Now we are ready to write down the equations of in-medium gluodynamics gen- 
eralizing Eq. ([3]) in the same way as Eqs. (fTUj) . (fTTI) are derived in electrodynamics. 
We introduce the nuclear permittivity and denote it also by e since it will not lead 
to any confusion. After that one should replace E a in Eqs. (JSJ), (J9]) by eE a and get: 

e(divE a - gf abc A b E c ) = p a , (13) 
dE 

curlB a - e-^- - gf abc (e<5> b E c + [A 6 B C ]) = j a . (14) 

The space-time dispersion of e is neglected here. 

In terms of potentials these equations are cast in the form: 

d 2 A 1 d 

A A a - = -Ja - ^/abc(-(curl[A b , A c ] + [A b curlA c ]) + 7^(A fe $ c ) - 

dA 1 
e ^^-Qf ~ e$bgrad$ c - -gf cmn [A b [A m A n }} + gef cmn $ b A m $ n ), (15) 



A $ a - 6^ = -^+ ^ a6c (2A 6 grad$ c + A b ^ - e^ c ) + 
at 1 e ot at 



g 2 famnfnlbAmM$b- (16) 
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If the terms with explicitly shown coupling constant g are omitted, one gets the 
set of Abelian equations which differ from electrodynamical equations fllOp . (II ip by 
the color index a only. Their solutions are shown in the next section. The external 
current is ascribed to a parton fast moving relative to other partons " at rest" . The 
crucial distinction between Eq. Q and Eqs. ffT5j) . ffTBT) is that there is no radiation 
(the field strength is zero in the forward light-cone and no gluons are produced) in 
the lowest order solution of Eq. ([3]) and it is admitted for Eqs. ffTol) . f[T5"j) because e 
takes into account the collective response (polarization) of the nuclear matter. 

The omitted above terms are of the order of g 3 because the potentials and the 
classical current J M are linear in g. They can be taken into account as a perturbation. 
It was done in [T31 [33] for in- vacuum gluodynamics. For in-medium gluodynamics 
they are considered in Section 7 and Ref. [15J. 

3 Cherenkov gluons as the classical lowest order 
solution of in-medium gluodynamics 

The classical solution of Eqs. (fT5|) . ({TBI) immediately leads to the notion of Cherenkov 
gluons at e > 1 in analogy with Cherenkov photons in electrodynamics. The unique 
feature is independence of the coherence of subsequent emissions by an external 
current on the time interval between these processes. 

The problem of the coherence length for Cherenkov radiation was extensively 
studied (161 El- It was shown that the tu-component of the field of a current can 
be imitated by a set of oscillators with frequency ui situated along the trajectory. 
The waves from all oscillators add up in the direction given by the Cherenkov angle 
9 independent on the length of the interval filled in by these oscillators. The phase 
disbalance A<p between emissions with frequency ui = k/y/e separated by the time 
interval At (or the length Az = vAt) is given by 

A0 = uAt - kAz cos 9 = kAz(-^-= - cos 9) (17) 

up to terms which vanish for large distances between oscillating sources and the 
detector. For Cherenkov effects the angle 9 is 

cos9=-^. (18) 
Vy/e 

The coherence condition A0 = is valid independent of Az. This is a crucial prop- 
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erty specific for Cherenkov radiation onljj^. Thus the change of color at emission 
vertices is not important if one considers a particular a-th component of color fields 
produced at Cherenkov angle. Therefore the fields ($ a , A a ) and the classical current 
for in-medium gluodynamics can be represented by the product of their electrody- 
namical expressions ($, A) and the color matrix T a . As a result, one can neglect the 
"rotation" of color at emission vertices and use in the lowest order for Cherenkov 
gluons the well known formulae for Cherenkov photons just replacing a by QsCa 
for gluon currents in probabilities of their emission. Surely, there is radiation at an- 
gles different from the Cherenkov angle (fl8l) . For such gluons one should take into 
account the coherence length and color rotation considering corresponding Wilson 
lines [IE]. 

Let us remind the explicit Abelian solution for the current with velocity v along 
z-axis 

j(r,i) = vp(r,t) = 4ngv6(r-vt). (19) 

In the lowest order the solutions for scalar and vector potentials are related so 
that 

A (1) (r,t) = ev$ (1) (r,t), (20) 

where the superscript (1) indicates the solutions of order g. 

Therefore the explicit expressions for $ suffice. Using the Fourier transform, the 
lowest order solution of Eq. ffTTT) with account of f[T9~j) can be cast in the form 

The integration over the angle in cylindrical coordinates gives the Bessel function 
Jo{k±r±). Integrating over the longitudinal component k z with account of the poles 
due to the denominatoi@ and then over the transverse one k±, one gets the following 
expression for the scalar potential 



d>«(r, t) = 29 ^-z-^V^T) 
e ^{vt~ z) 2 -r\{ev 2 -I) 

Here r± = ' x 2 + y 2 is the cylindrical coordinate, z is the symmetry axis. The 
cone 

z = vt- r ± Vev 2 - 1 (23) 



3 The requirement for Atfi to be a multiple of 2ir (or a weaker condition of being less or of 
the order of 1) in cases when Cherenkov condition is not satisfied imposes limits on the effective 
radiation length as it happens, e.g., for Landau-Pomeranchuk or Ter-Mikaelyan effects. 

4 These poles are at work only for Cherenkov radiation! 
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determines the position of the shock wave due to the ^-function in Eq. (1221) . The 
field is localized within this cone. The Descartes components of the Poynting vector 
are related according to Eqs. (1221) . (120]) by the formulae 

(z - vt)x q q (g - vt)y , , 

'_!_ '_L 

so that the direction of emitted gluons is perpendicular to the cone ( |23i) and defined 
by the Cherenkov angle 

tan 2 fl = = ev * _ i, (25) 

which coincides with (I18p . 

The higher order terms (g 3 ...) can be calculated using Eqs. (JT5|) . (TTB1) (see 
Section 7). 

The expression for the intensity of the radiation is given by the Tamm-Frank 
formula (up to Casimir operators) 

dW f 1 

-^-=4vra 5 / udw 1-— ). (26) 
dz J v 2 e 

It is well known that it leads to infinity for constant e. The (^-dependence of e (its 
dispersion) usually solves the problem. 

For absorbing media e acquires the imaginary part e = 6\ + ie^- The sharp front 
edge of the shock wave is smoothed. The angular distribution of Cherenkov radiation 
widens. The ^-function at the angle (TTSi) . fTSBT) is replaced by the a'la Breit-Wigner 
angular shape [2U] with maximum at the angle given by (1181) . (1251) but with |e| in 
place of e and the width proportional to the imaginary part. 

The energy loss dW per the length dz is determined by the formula 



dz 

In the lowest order 

J4 



-gE z . (27) 



EP = i I ^M«(k, W ) - k z ^(k,u)]e^\ (28) 
and in ^-representation one gets in QCD 

u 2 e(ev % — 1) 
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4g = (30) 
C = cos#. (31) 

The subscript of e a is omitted because there will be no summation on it even at 
higher order (see Section 7). 

Thus, for ultrarelativistic radiating centers, the differential Cherenkov energy 
loss spectrum in the opaque medium per length dz reads [201 [211 122] 

1 dW _ 4a s C v cosfl(l-cos 2 g)lV) 
udz du d<f)dcos0 ~ tt (cos 2 6 - CM) 2 + T 2 (u) ' 

where 0, are polar and azimuthal angles with respect to the emitter propagation 
direction, 



CM 



(33) 



so that C( w ) controls the location of the maximum and T(u) controls the spreading 
around it. For (e 2 /e!) 2 <C 1, the energy-dependent maximum of the differential 
spectrum f l32|) is at 

cos 2 6» max (u;) « „ v— 27 r - (34) 

Without absorption, the potential (l22il is infinite on the cone. With absorption, it 
is finite everywhere except the cone vertex and is inverse proportional to the distance 
from the vertex. Absorption induces also longitudinal excitations (chromoplasmons) 
which are proportional to the imaginary part of e and usually small compared to 
transverse excitations. 

Several experimental observations may be explained as stemming from emission 
of Cherenkov gluons. Most accurate are the RHIC data with high statistics about the 
two-humped structure of particle distributions around the away-side jets in central 
nucleus- nucleus collisions at y/s = 200 GeV per nucleon [H [21 [31 HI [6] . Their fit 
allowed to get values of the real and imaginary parts of the nuclear permittivity 
[22J. The earlier cosmic ray event [23], which initiated the idea about Cherenkov 
gluons, showed the smaller values of the nuclear permittivity [241 [25]. Quite small 
values of e were also obtained from low-statistics samples of nuclei interactions at 
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SPS energies [23 EH]- This indicates on energy dependence (dispersion) of e. 
The SPS data |29j on asymmetry of the shapes of resonances passing through the 
nuclear medium can also be interpreted in terms of Cherenkov gluons [30J. Let us 
start with fits of the RHIC data. 



In Eq. f l32|) the angle 6 of emission of Cherenkov gluons is measured with respect to 
the direction of propagation of the radiating color current. In the experimental setup 
used at RHIC the special trigger method has been used. Among all central Au-Au 
collisions at 200 GeV per nucleon there were chosen only those where the partons 
of colliding nuclei scattered at the angle about rr/2. One of them passed a very thin 
layer of a nucleus and produced the jet similar to those in pp-collisions. Another 
one traversed through the whole nucleus and formed the away-side jet. Namely 
around this jet the rings of hadrons were detected and interpreted as originating 
from Cherenkov gluons. The partons initiating the hard away-side jets play the role 
of the radiating color currents. Since the trigger is placed at 7r/2 to the colision axis 
of initial partons, the trigger parton is detected at this angle. The away-side parton 
is also at this angle in the opposite direction if the energies of colliding partons 
are equal. It was shown in [25] that the mismatch of these energies is unimportant 
because the structure functions decrease fast for larger mismatch. The background 
from gluons radiated by the aside moving partons is low and can result in a slight 
widening of observed humps. Due to the forward-backward symmetry it does not 
influence the positions of the maxima. Thus, in the first approximation, for analyt- 
ical estimates one can consider only those partons whose direction of propagation is 
orthogonal to the collision axis z. 

In order to compare the spectrum (1321) with the experimental data one should 
rewrite it in terms of experimental laboratory polar and azimuthal angles 6l and 
<Pl- It is easy to see that for the above geometry 



for 4>l counted from the away-side jet. 

The number dN of emitted gluons per length dz with energy and angles within 
intervals du dO^ dcf)^ is 



4 Comparison with RHIC data 



cos# 



sin 6l cos 4>l\, 



dN 



4asCy | sinOi cos0l|(1 — sin 2 9l cos 2 ^>l)Y{u) 



(35) 



dz dui d(j)L d cos 9 l 



(sm 2 e L cos 2 L -cH) 2 + r 2 {u) 
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The azimuthal gluon distribution dN/ difii is obtained from Eq. 035)) by integrat- 
ing over 9 1. This can be done analytically: 




where 

A(u,<j> L ) = r 2 (w)+C 2 M-CMcos 2 L , B(u,<j> L ) = I»cos 2 L . (37) 

The gluon spectrum (|36|) reveals the double-humped structure. The original flow 
of Cherenkov glue described by Eq. (l3l)j) should, of course, transform into that of 
final hadrons. This can be done only with the help of Monte-Carlo models. 

To describe properly the angular pattern observed in the flow of final hadrons in 
high energy nuclear collisions [TJ [6j [31 H] in terms of the Cherenkov gluon radiation 
one has to 

(a) Describe a kinematical pattern characterizing the initially produced hard 
partons serving as colored sources of gluon Cherenkov radiation taking into account 
the experimental cuts on pseudorapidity and transverse momentum [6JH]. 

(b) Write down a spectrum of Cherenkov gluon radiation in the opaque medium 
for the above-described color sources. 

(c) Describe a conversion of Cherenkov gluons into observed hadrons taking into 
account the experimental cuts on transverse momenta of final hadrons [21 H] • 

These three steps are realized in a Monte-Carlo procedure. In total the model 
includes three parameters, described in more details below - two related to the 
Cherenkov gluon radiation, the magnitudes of the real e\ and imaginary 62 parts of 
the permittivity of the medium, and one related to rescattering of Cherenkov gluons 
in the process of hadronization Aj_. The values of these parameters are fitted to 
achieve the best possible agreement with the experimental spectra [6JH]. 

The description of the initial hard color sources was performed through a Monte 
Carlo simulation of hard parton-parton scattering at RHIC energies with PYTHIA 
|31j . This authomatically allowed to take into account possible mismatch of energies 
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of initial colliding partons. An initial pool of two-parton configurations resulting 
from these scatterings consisted of those in which the trigger jet gave rise to a 
pion hitting the trigger intervals in rapidity, | r]\ll aT — 0.7 and | ?7|p r henix < 0.35, and 
transverse momentum, 3 GeV < | p^ b | < 4GeV for STAR and 2 GeV < |pj*| < 
3 GeV for PHENIX [51 H] . The fragmentation of gluon (quark) generating the trigger 
near-side jet into pions was described by standard fragmentation functions [32J. 

The away-side jet gave rise to a pion hitting the rapidity intervals | r]\^ ar < 1 
and | rj\ ph enix < 0.35, and transverse momentum, 1 GeV < | p^j < 2 GeV for STAR 
and 2 GeV < I p^ b | < 3 GeV for PHENIX [5J Hj. For simplicity we have considered 
only the dominating subset of events with gluonic away-side jet. 

The angular distribution of emitted Cherenkov gluons is given by the expresion 
( I3"5|) . Let us stress once again that in our Monte-Carlo procedure we consider a 
more general situation where the direction of the gluon generating the away-side jet 
is fixed within each configuration satisfying the above-described trigger conditions 
imposed on the properties of the near-side jet in experiment. 

To realize a Monte-Carlo procedure of generating the Cherenkov spectrum we 
have to specify the functions e\{uj) and €2{uj). Taking into account that the experi- 
mental cuts on (laboratory frame) intervals of transverse momenta of the away-side 
hadrons are very strict, 5p^ b | < 1 GeV both for STAR [5J and PHENIX [I] we 
can with a good accuracy neglect the effects of dispersion and consider energy- 
independent ei 2 = const. The values of these constants are determined through 
fitting the experimental data. Within this assumption the spectrum of produced 
Cherenkov glue is simply energy- independent, dN/dudl = const. Realistically the 
possibility of Cherenkov emission is restricted to some finite interval of energies 
uj < c<j max so that 



where kv ax is the highest characteristic resonance excitation energy of the medium_|. 
In our simulations we have chosen u; max = 3.5 GeV in accordance with the upper 
limits of experimental intervals. We have verified that our results are in fact weakly 
sensitive to the exact value of u; max . 

The thus generated flow of Cherenkov glue should, of course, be transformed into 
that of final hadrons. There exist several phenomenological schemes describing this 
conversion. In our case it is convenient to use a language of fragmentation functions 
Dg(x, p± | Q 2 ) which generically describe a probability for a gluon with energy E and 

5 For a discussion of interrelation between Cherenkov gluons and hadronic resonances see [30j 
and section 6. 
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virtuality Q 2 to produce a hadron with the energy xE and transverse momentum 
p± measured with respect to the direction of propagation of the initial gluon. 

One is first tempted to ascribe the same shape to the spectrum of Cherenkov 
hadrons relying on the soft blanching hypothesis, equivalent to assuming Dg(x, p_i_ | Q 2 ) oc 
5(1 — x) supported by the experimental evidence obtained in the physics of multi- 
particle production from e + e~ — ► jets. However, this does not lead to the fully 
satisfactory description of experimental data because of predicting, in contradiction 
with the experimental data, that the probability of radiating Cherenkov glue at 
angles \4>l\ > 7r /2 with respect to the away-side jet is strictly zero. These "dead 
zones" appear because the Cherenkov gluon "halo" can not spread to angles larger 
than 7r/2 to the direction of the emitter. It is clear that such "dead zones" will in 
fact be present for any fragmentation mechanism that does not generate transverse 
momentum with respect to the direction of propagation of the initial Cherenkov 
gluon which in our case corresponds to considering a p_i_-independent fragmenta- 
tion function D^(x, \ Q 2 ). We confine our consideration to fragmentation to light 
hadrons because experimentally the share of other bosonic resonances is quite small 
(in particular, p : to : 0=10:1:2 according to [2S])- In our analysis we used a simple 
parametrization of the fragmentation function Dg(x, p±\ Q 2 ): 

^(x, P jg 2 )cx(l-x) 3 - 7 =i==exp|-^} , (39) 

where the x-dependence is that of a gluon fragmentation function at the reference 
scale Qo = 2 GeV [32]. Phenomenologically, the parameter Aj_ accounts for trans- 
verse momenta acquired by pions both due to the decay of intermediate resonances 
and because of rescattering processes. 

Thus, we have a set of three parameters: ei, e 2 and a fraction Aj_. 

The values of the parameters that provide the best description of experimental 
data of STAR and PHENIX collaborations are shown in Table 1. 



Table 1 



Experiment 


^max 


ei 


^2 


A ± 


STAR 


1.04 rad 


5.4 


0.7 


0.7 GeV 


PHENIX 


1.27 rad 


9.0 


2.0 


1.1 GeV 



Let us stress that a difference in the fitted values of e\ for STAR and PHENIX 
originates from different positions of angular maxima ^ max in the corresponding 
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experimental data. (May be, it could indicate on necessity to take the dispersion 
into account.) However these values are quite stable in both above approaches 
because they are defined by maxima positions but not by the humps widths. 

The value of e 2 , on the contrary, is influenced by the widths. Nevertheless, the 
ratio (e 2 /ei) 2 < 0.05 <C 1 stays quite small in all cases. 

The transverse smearing parametrized by Aj_ looks very reasonable at the hadronic 
scale. 

The resulting angular spectra for STAR and PHENIX are shown in Figs. 1 
and 2 correspondingly. We see that the positions of the maxima (and therefore the 
values of ei) are quite stable to accounting for additional smearing on top of that 
in Eq. (I36p . However, it is indeed important in achieving a good description of 
the widths of humps in experimental data as seen in Figs. 1 and 2. The shape of 
humps in the former "dead zone" determines the parameter Aj_, which, in its turn, 
influences e 2 . 

To conclude, the values of the real part of the nuclear permittivity e\ found 
from the fit to experimental data of RHIC (albeit somewhat different in the two 
experimental sets) are determined with good precision. Their common feature is 
that they are noticeably larger than 1. This shows that the density of scattering 
centers is quite large and the nuclear medium reminds a fluid rather than a gas (for 
more details see [25]). The accuracy of the estimate of the values of e 2 is much less 
but more important is the fact that they are rather small compared with e\. 

5 Cosmic ray events and SPS data 

The RHIC experiments with high statistics used a special geometry of events and 
corresponding triggers to detect rings of hadrons around the away-side jets. Lower 
statistics data obtained at SPS and in cosmic rays also gave some indications on 
ring-like events. The very first events with the ring-like structure were observed in 
non-trigger cosmic ray experiments. Namely they initiated the idea about Cherenkov 
gluons [8]. Two rings more densely populated by particles than their surroundings 
were noticed in the cosmic ray event [23J initiated by a primary with energy about 
10 16 eV close to LHC energies. It is demonstrated in Fig. 3 where the number of 
produced particles is plotted as a function of the logarithm of the distance from 
the collision axis (proportional to pseudorapidity at small angles). It clearly shows 
two maxima. This event has been registered in the detector with nuclear and X-ray 
emulsions during the balloon flight at the altitude about 30 km. Approximately at 
the same time the similar event with two peaks was observed at 10 13 eV [31]. Two 
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colliding nuclei must give rise to two peaks. Some events with one peak (due to the 
limited acceptance of the installation?) were shown even earlier [35, 36J. The peaks 
were interpreted as effects due to Cherenkov gluons emitted by the forward and 
backward moving initial high energy partons. When the two-dimensional distribu- 
tion of particles was considered in the azimuthal plane (called as target diagram in 
cosmic rays experiments), this event revealed two (forward and backward in c.m.s.) 
densely populated ring-like regions within two narrow intervals of polar angles (cor- 
responding to peaks in Fig. 3) but widely distributed in azimuthal angles for each 
of them. Therefore such events were called as ring-like events. 

It was shown that the ring-like structure can be revealed by the event-by-event 
wavelet analysis [371 EH] • 

The argument in favor of high energy Cherenkov gluons stems from the behavior 
of the real parts of hadronic forward scattering amplitudes KeF(E, 0°). The intensity 
of the Cherenkov effect is proportional to the excess of the permittivity e or of the 
refractive index n over 1. They are related by the formula e = n 2 . The general 
relation of the scattering theory (see, e.g., [39] ) between the refractive index and the 
forward scattering amplitude F(E, 0°) states that 



Here E denotes the energy, N is the density of the scatterers (inhomogeneities) of the 
medium, is the number of scatterers within a single nucleon, m n the pion mass, 
<y(E) the cross section and p(E) the ratio of real to imaginary parts of the forward 
scattering amplitude F(E). Thus the emission of Cherenkov gluons is possible only 
for processes with positive KeF(E) or p(E). As follows from (HOI) , the values of 
An can be completely different at low and high energies because of different values 
of N and KeF(E, 0°) involved. Unfortunately, we are unable to calculate directly 
in QCD KeF(E, 0°) for gluons and have to rely on analogies and our knowledge 
of properties of hadrons. The only experimental facts we get about this medium 
are brought by particles registered at the final stage. They have some features in 
common which (one can hope!) are also relevant for gluons as the carriers of the 
strong forces. These are the resonant behavior of hadronic amplitudes at rather low 
energies and positiveness of KeF(E, 0°) for all measured hadronic processes at very 
high energies. Thus we wait for Cherenkov effects in the resonance region and at 
high energies. 

The real parts of the resonance amplitudes are positive in their low-mass wing. 
This fact will be extensively used in the next section. Here, it is important that the 
energies of particles produced within the humps in RHIC data are rather low on 



An = Ken — 1 



2nNReF(E, 0°) _ 3m> 
E 2 ~ 2E 2 



ReF(E) 



3m^z/ 
8nE 



a(E)p(E). (40) 
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the scale of initial energies. Even the energies of jets-parents are lower than 5 GeV. 
Therefore the whole effect at RHIC can be related to the low-energy region. 

In cosmic ray events the energies of partons are high and the effect can be 
ascribed to the high-energy behavior of KeF(E, 0°). 

One of the most intriguing problems is to understand properly the fact that the 
RHIC and cosmic ray data were fitted with very different values of the refractive 
index close to 3 and 1, correspondingly. This could be interpreted as due to the 
difference in values of x (the parton share of energy) and Q 2 (the transverse mo- 
menta). It is well known that the region of large x and Q 2 corresponds to the dilute 
partonic system. At low x and Q 2 the density of partons is much higher. 

The density of scattering centers at RHIC conditions is very high as estimated 
from Eq. (I4(jp [25J, about tens partons per proton volume. Here one deals with 
rather low x and Q 2 . Therefore, one should expect the large density of partons in 
this region and high An ~ 1. It is interesting to note that the two-hump structure 
disappears in RHIC data at higher p t where the parton density must get lower. It 
corresponds to smaller n and 8, i.e. humps merge in the main away-side peak. 

In the cosmic ray event one observes effect due to leading partons with large x. 
Also, the experimentalists pointed out that the transverse momenta in this event 
are quite large. In this region one would expect for low parton density and small 
An < I. 

The ring-like events with values of e close to but larger than 1 were also found 
from the SPS non-trigger experiments [27J, 12S]- 

Thus the same medium can be probably seen as a liquid or a gas depending on 
the parton energy and transferred momenta. This statement can be experimentally 
verified by using triggers positioned at different angles to the collision axis and 
considering different transverse momenta. In that way, the hadronic Cherenkov 
effect can be used as a tool to scan (1/x, Q 2 )-plane and plot on it the parton densities 
corresponding to its different regions. We discuss this problem also in Section 8. 

6 Asymmetry of shapes of resonances produced 
inside the nuclear medium 

The low-energy Cherenkov gluons may be at the origin [30J of another interesting 
effect - the asymmetry of shapes of resonances created inside the nuclear medium. 
There exist numerous experimental data [29l S21 S3, SSJ, STJ about 

the in-medium modification of widths and positions of prominent vector-meson res- 
onances. They are mostly obtained from the shapes of dilepton mass and transverse 
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momentum spectra in nucleus-nucleus collisions. Such in-medium effects were tied 
theoretically to chiral symmetry restoration a long time ago [48J. 

A significant excess of low-mass dilepton pairs yield for p meson over expectations 
from hadronic decays is observed in the high-statistics experiment [29] . 

Several approaches have been advocated for explanation of the excess. Strong 
dependence of the parameters of the effective Lagrangian on the temperature and 
the chemical potential was assumed in [49"| I5U]. The hydrodynamical evolution was 
incorporated in [51] to describe the spectra. The QCD sum rules and dispersion 
relations have been used [52j [53] to show that condensates decrease in the medium 
leads to both broadening and slight downward mass shift of resonances. The similar 
conclusions have been obtained from more traditional attempts using either the 
empirical scattering amplitudes with parton-hadron duality [541 [55] or the hadronic 
many-body theory [561 [571 EH] • 

In the latest approach, which pretends to provide the best description of exper- 
imental plots, the in-medium V-meson spectral functions are evaluated. The excess 
of dilepton pairs below p-mass is ascribed to antibaryonic effects. This conclusion is 
the alternative to more common ideas about the chiral restoration at high energies. 
It asks for some empirical constraints to fit the observed excess. 

We proposed [30] another possible source of low-mass lepton pairs. Namely, the 
emission of Cherenkov gluons may provide a substantial contribution to the low 
mass region for any resonance. 

Qualitatively, the observed low mass excess of lepton pairs is easy to ascribe to the 
gluonic Cherenkov effect if one reminds that the index of refraction of any medium 
exceeds 1 within the lower wing of any resonance (the p-meson, in particular). 

This feature is well known in electrodynamics (see, e.g., Fig. 31-5 in [59]) where 
the atoms behaving as oscillators emit as Breit-Wigner resonances when get ex- 
cited. This results in the indices of refraction larger than 1 within their low-energy 
wings. In QCD, one can imagine that the nuclear index of refraction for gluons 
in the hadronic medium behaves similarly in the resonance regions. In classical 
electrodynamics, it is the dipole excitation of atoms in the medium by light which 
results in the Breit-Wigner shape of the amplitude F(E,0°). In hadronic medium, 
there should be some modes (quarks, gluons or their preconfined bound states, con- 
densates, blobs of hot matter...?) which can get excited by the impinging parton, 
radiate coherently if n > 1 and hadronize at the final stage as hadronic resonances 

mii]. 

The scenario, we have in mind, is as follows. Any parton, belonging to a colliding 
nucleus, can emit a gluon. On its way through the nuclear medium the parton 
excites some internal modes. Therefore it affects the medium as an "effective" wave 
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which accounts also for the waves emitted by other scattering centers (see, e.g., 
|39j). Beside incoherent scattering, there are processes which can be described as 
the refraction of the initial wave along the path of the coherent wave. The Cherenkov 
effect is the induced coherent radiation by a set of scattering centers placed on the 
way of propagation of a gluon. That is why the forward scattering amplitude plays 
such a crucial role in formation of the index of refraction. At low energies its excess 
over 1 is related to the resonance peaks as dictated by the Breit-Wigner shapes 
of the amplitudes (for the similar well known explanation in electrodynamics see, 
e.g., [59]). In experiment, usual resonances are formed during parton interactions 
and subsequent color neutralization process. Up to now there is no quantitative 
explanation of resonant amplitudes at the level of QCD-partons. One can only admit 
some attractive forces acting within definite energy intervals which lead to resonant 
behavior. In addition to these forces between the individual partons the collective 
effects can appear in the nuclear medium. As discussed above, there are several 
indications on such effects in RHIC experiments like J/ ^-suppression, azimuthal 
asymmetry t> 2 , jet quenching, Cherenkov rings etc (see, e.g., [H [21 12 E]). Among 
them, Cherenkov gluons are distinguished as the result of the fully coherent reaction 
of the medium. Only in this case the phase of the amplitude does not depend on the 
time interval between emissions (see Section 4). The necessary condition n > 1 is 
satisfied within the left-wing resonance region and coherent effect at low energies can 
be observed only there. Similarly to other partons, Cherenkov gluons participate in 
formation of resonances contributing to their low-mass wing (and, in particular, to 
dilepton spectra at these energies). Therefore these radiative effects can add to the 
ordinary ones (strings?) namely in this specific region. This contribution should 
be proportional to An according to the theory [16J. The shock wave formed by 
Cherenkov gluons [7] pushes out the hadronic states just with masses within the 
low-mass wings of resonances. Let us note that the common string description of 
the processes does not take into account collective effects. 

Thus, the ordinary Breit-Wigner shape of the cross section for resonance pro- 
duction must be modified by the coherent medium response. The easiest way to 
observe this effect would be by measuring the dilepton mass spectra for resonances. 
For all of them, they would acquire the additional term proportional to An at masses 
below the resonance peak. According to Eq. (f4"U|) it must be proportional to the 
real part of the Breit-Wigner amplitude with the ratio of real to imaginary parts of 
Breit-Wigner amplitudes equal to 




(41) 
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Here M is the total c.m.s. energy of two colliding objects (the dilepton mass), 
m p =775 MeV is the in-vacuum p-meson mass. Thus its M-dependence is well 
defined. This term vanishes for M > m p because only positive An lead to the 
Cherenkov effect. Namely it describes the distribution of masses of Cherenkov states. 
Therefore, the shape of the mass distribution near the p-meson can be described 



The first term corresponds to the ordinary Breit-Wigner cross section with a constant 
normalization factor A. According to the optical theorem it is proportional to the 
imaginary part of the forward scattering amplitude. The second term is due to the 
coherent response of the medium and its weight relative to ordinary processes^ is 
described by the only adjustable parameter w. As for any classical effect, one can 
not calculate the cross section for Cherenkov effect and determine the parameter w 
theoretically. Therefore, we leave w as an adjustable parameter. The only hope is 
that it is not very small because such an effect has been observed recently in high 
energy heavy-ion collisions at RHIC (section 4) as the nimbus around away-side jets 
with rather high probability. Our fit below supports the assumption that the value 
of w is quite noticeable and corresponds to higher energy observations. 

In these formulas, one should take into account the in-medium modification of 
the height of the peak and its width. In principle, one could consider free 
in-medium parameter as well. We rely on experimental findings that its shift in the 
medium is small. All this asks for some dynamics to be known. In our approach, 
it is not determined. Therefore, first of all, we just fit the parameters A and T by 
describing the shape of the mass spectrum at 0.75 < M < 0.9 GeV measured in [29J 
and shown in Fig. 4. In this way we avoid any strong influence of the 0-meson. Let 
us note that w is not used in this procedure. The values A=104 GeV 3 and T = 0.354 
GeV were obtained. The width of the in-medium peak is larger than the in-vacuum 
p-meson width equal to 150 MeV. It agrees with previous findings of widening of 
in-medium p-mesons. 

Thus the low mass spectrum at M < m p depends only on a single parameter 
w which is determined by the (theoretically unknown) relative role of Cherenkov 
effects and ordinary mechanism of resonance production. It is clearly seen from Eq. 

6 We consider only p-mesons here. To include other mesons, one should evaluate the similar 
expressions. 

7 We assume it to be independent of M because the relative probability of Cherenkov emission 
to ordinary processes can hardly change within this rather short interval of masses. 



by the following formula[§ 




(42) 
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( }4"2"j) that the role of the second term in the brackets increases for smaller masses M. 
The excess spectrum in the mass region from 0.4 GeV to 0.75 GeV has been fitted by 
w = 0.19. The slight downward shift about 40 MeV of the peak of the distribution 
compared with m p may be estimated from Eq. (|42p at these values of the parameters. 
This agrees with the above statement about small shift compared to m p . The total 
mass spectrum (the dashed line) and its widened Breit-Wigner component (the solid 
line) according to Eq. (T4"2l with the chosen parameters are shown in Fig. 4. The 
overall description of experimental points seems quite satisfactory. The contribution 
of Cherenkov gluons (the excess of the dashed line over the solid one) constitutes 
the noticeable part at low masses. The formula fj4*2l) must be valid in the vicinity of 
the resonance peak. Thus we use it for masses larger than 0.4 GeV only. 

From general principles one would expect slightly lower for low-mass dilepton 
pairs from coherent Cherenkov processes than for incoherent scattering. Qualita- 
tively, this conclusion is supported by experiment [29]. 

Whether the in-medium Cherenkov gluonic effect is as strong as shown in Fig. 
4 can be verified by measuring the angular distribution of the lepton pairs with 
different masses. The trigger-jet experiments similar to that at RHIC are necessary 
to check this prediction. One should measure the angles between the companion jet 
axis and the total momentum of the lepton pair. The Cherenkov pairs with masses 
between 0.4 GeV and 0.7 GeV should tend to fill in the rings around the jet axis. 
The angular radius 9 of the ring is determined by the usual condition. 

Another way to demonstrate it is to measure the average mass of lepton pairs 
as a function of their polar emission angle (pseudorapidity) with the companion jet 
direction chosen as z-axis. Some excess of low-mass pairs may be observed within 
the rings. Baryon-antibaryon effects can not possess signatures similar to these ones. 

Some indications on the substructure with maxima at definite angles have been 
found at the same energies by CERES collaboration [52]. It is not clear yet if it 
can be ascribed to Cherenkov gluons. To recover a definite maximum, it would be 
better to detect a single parton jet, i.e. to have a trigger. 

The prediction of asymmetrical in-medium widening of any resonance at its 
low-mass side due to Cherenkov gluons is universal. This universality is definitely 
supported by experiment. Very clear signals of the excess on the low-mass sides of 
p, uj and mesons have been seen in (431 E] • This effect for uj- meson is also studied 
in [45]. Slight asymmetry of 0-meson near 0.9 - 1 GeV is noticeable in the Fig. 
4 shown above but the error bars are large there. We did not try to fit it just to 
deal with as small number of parameters as possible. There are some indications at 
RHIC (see Fig. 6 in [45]) on this effect for J/ip-meson. However, the accuracy of 
RHIC data both for p and J / ip is not enough to get quantitative conclusions. 
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To conclude, the new mechanism is proposed for explanation of the low-mass 
excess of dilepton pairs observed in experiment. It is the Cherenkov gluon radiation 
which adds to the ordinary processes at the left wing of any resonance. Only one 
adjustable parameter w is used to describe its contribution to the dilepton spectra. 
The universal nature of their asymmetry for all resonances is predicted. 

7 Higher order effects leading to the "color rain- 
bow" 

It is interesting to note that the general first order solutions of Eqs. (1151) . (fl6l) . 
explicitly shown by Eqs. (I22I) . (1201) for Cherenkov gluons, stay valid at all orders 
if the nuclear permittivity does not depend on color. The higher order terms con- 
tribute only if the permittivity differs for different colors. The general procedure 
of their calculation is the same as for in- vacuum QCD jTOj, [131 EH]- After getting 
explicit lowest order solution one exploits it together with the non-Abelian current 
conservation condition to find the current component proportional to g 3 . Then with 
the help of Eqs. f[T5"j) . ffTB]) one finds the potentials up to the order g 3 . They can be 
represented as integrals convoluting the current with the corresponding in-medium 
Green function. The even higher order corrections may be obtained in the same 
way. We outline briefly the path to the next order effects and qualitative results. 

The third order terms of the potentials can be found [15] from the following 
expressions 

e(A*?> - e^) = -gf abc ie}(e f v* - 1) + (2c + e>> 6c ] , (43) 

A 6c = (e b - e c )/2; e f = (e b + e c )/2, 

f) 2 A {3) (9<J> (1) 
A Afl - = -gf abc [e f (e f v - l/v) + A bc (2v(e + e f ) - l/„)]-^-«|>« (44) 



The density of the energy loss is proportional according to Eq. (1271) to 

Efl = < / (^i (W 3 l(k,c) - k(^\k,co)] - igf abc J ^ A W{k>)*V{k - fc')) 

(45) 

The last term corresponds to the truly non-abelian correction. Only the linear 
terms in A bc must be left everywhere because of the antisymmetry of f abc . It demon- 
strates that the classical first order solution stays valid for the color independent 
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permittivity. At c can differ from either due to dispersion effects or because of 
explicit dependence of permittivity on color. 

The calculation of the spectra proceeds [15] in the same way as it was done above 
for first order terms. For the sake of brevity, we reproduce here only the first term 
in the limit e — 6f, v — 1: 



dN [A) 

an a{A<$>)lim 

dldxdu 



3g| 
4v/i 



fabcQbQclKL 



1) 1/2 A 



be 



e 2 (ex 



!)3/2 



(46) 



The terms oc 



and oc 



appear which are differ- 



[(x-xo^ + iuxo) 2 ] 5 /* allyx ^ y/x-[(x-x ) 2 + {vxa) 2 ? /A 

ent (albeit somewhat similar) to the lowest order terms. Now, the shape is different 
from simple a'la Breit-Wigner form. The typical proportionality to An is seen. 

To conclude, one can say that in the case of Af, c ^ the non-abelian color 
quantum rainbow appears due to the higher order terms. 



8 The nuclear permittivity and the rest system of 
the nuclear matter 

In electrodynamics the permittivity of real substances depends on uj. Moreover it 
has the imaginary part determining the absorption. E.g., Ree for water (see [63J) 
is approximately constant in the visible light region (y/e ~ 1-34), increases at low 
uj and becomes smaller than 1 at high energies tending to 1 asymptotically. The 
absorption (Im e) is very small for visible light but dramatically increases in nearby 
regions both at low and high frequencies. Theoretically this behavior is ascribed to 
various collective excitations in the water relevant to its response to radiation with 
different frequencies. Among them the resonance excitations are quite prominent 
(see, e.g., [59]). Moreover, the medium considered is at rest and the permittivity 
values are determined just in this frame. Even in electrodynamics, the quantitative 
theory of their energy dependence is still lacking, however. 

Then, what can we say about the nuclear permittivity and the frame to define 
it in? 

The partons constituting high energy hadrons or nuclei interact during the colli- 
sion for a very short time. Nevertheless, there are experimental indications that an 
intermediate state of matter (CGC, QGP, nuclear fluid ...) is formed and evolves. 
Those are J/^-suppression, jet quenching, collective flow (t> 2 ), Cherenkov rings of 
hadrons etc. They show that there is collective response of the nuclear matter to 
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color currents moving in it. Unfortunately, our knowledge of its internal excitation 
modes is very scarce, much smaller than in electrodynamics. 

The permittivity is the internal property of a medium which demonstrates the 
medium response to the induced current. Its quantitative description poses prob- 
lems even in QED. It becomes more difficult task in QCD where confinement is not 
understood. The attempts to calculate the nuclear permittivity from first principles 
are not very convincing. It can be obtained from the polarization operator. The 
corresponding dispersion branches have been computed in the lowest order pertur- 
bation theory [Ml l65l [66] . Then the properties of collective excitations have been 
studied in the framework of the thermal field theories (for review see, e.g., |67j). 
Their results with additional phenomenological ad hoc assumption about the role 
of resonances were used in a simplified model of scalar fields [H] to show that the 
nuclear permittivity can be larger than 1 that admits Cherenkov gluons. 

Our main goal is to study the medium response to the external color current. 
We prefer to use the general formulae of the scattering theory [39] to estimate the 
nuclear permittivity. We have to rely on analogies and our knowledge of properties 
of hadrons. The only experimental facts we get about this medium are brought by 
particles registered at the final stage. They have some features in common which 
(one can hope!) are also relevant for gluons as the carriers of the strong forces. Those 
are the resonant behavior of amplitudes at rather low energies and positive real part 
of the forward scattering amplitudes at very high energies for hadron-hadron and 
photon-hadron processes as measured from the interference of the Coulomb and 
hadronic parts of the amplitudes. KeF(0°,E) is always positive (i.e., n > 1) within 
the low-mass wings of the Breit-Wigner resonances. This shows that the necessary 
condition for Cherenkov effects n > 1 is satisfied at least within these two energy 
intervals. This fact was used to describe experimental observations at SPS, RHIC 
and cosmic ray energies. The asymmetry of the p-meson shape at SPS [29] and 
azimuthal correlations of in-medium jets at RHIC 0, [3] were explained above by 
emission of comparatively low-energy Cherenkov gluons [30l[2lj . The parton density 
and intensity of the radiation were estimated. In its turn, cosmic ray data [23] at 
energies corresponding to LHC ask for very high energy gluons to be emitted by 
the ultrarelativistic partons moving along the collision axis [8J. Let us note the 
important difference from electrodynamics where n < 1 at high frequencies. For 
QGP the high-energy condition n > 1 is a consequence of its instability. 

The in-medium equations are not Lorentz- invariant. There is no problem in 
macroscopic electrodynamics because the rest system of the macroscopic matter is 
well defined and its permittivity is considered there. For collisions of two nuclei (or 
hadrons) it asks for special discussion. 



22 



Let us consider a particular parton which radiates in the nuclear matter. It 
would "feel" the surrounding medium at rest if momenta of all other partons (or 
constituents of the matter), with which this parton can interact, sum to zero. In 
RHIC experiments the triggers which registered the jets (created by partons) were 
positioned at 90° to the collision axis. Such partons should be produced by two 
initial forward-backward moving partons scattered at 90°. The total momentum 
of other partons (medium spectators) is balanced because for such geometry the 
partons from both nuclei play a role of spectators forming the medium. Thus the 
center of mass system is the proper one to consider the nuclear matter at rest in 
this experiment. The permittivity must be defined there. The Cherenkov rings 
consisting of hadrons have been registered around the away-side jet which traversed 
the nuclear medium. This geometry requires however high statistics because the 
rare process of scattering of initial partons at 90° has been chosen. 

The forward (backward) moving partons are much more numerous and have 
higher energies. However, one can not treat the radiation of such a primary parton 
in c.m.s. in the similar way because the momentum of spectators is different from 
zero i.e. the matter is not at rest. Now the spectators (the medium) are formed from 
the partons of another (target) nucleus only. Then the rest system of the medium 
coincides with the rest system of that nucleus and the permittivity should refer to 
this system. The Cherenkov radiation of such highly energetic partons must be 
considered there. That is what was done for interpretation of the cosmic ray event 
in [8]. This discussion clearly shows that one must carefully define the rest system 
for other geometries of the experiment with triggers positioned at different angles. 

Thus our conclusion is that the definition of e depends on the experiment geome- 
try. Its corollary is that partons moving in different directions with different energies 
can "feel" different states of matter in the same collision of two nuclei because of 
the dispersive dependence of the permittivity. The transversely scattered partons 
with comparatively low energies can analyze the matter with rather large permit- 
tivity corresponding to the resonance region while the forward moving partons with 
high energies would "observe" low permittivity in the same collision. This peculiar 
feature can help scan the (lnx, Q 2 )-plane as it is discussed in Section 5 and [25]. It 
explains also the different values of e needed for description of RHIC and cosmic ray 
data. 

These conclusions can be checked at LHC because both RHIC and cosmic ray 
geometry will become available there. The energy of the forward moving partons 
would exceed the thresholds above which n > 1. Then both types of experiments can 
be done, i.e. the 90°-trigger and non-trigger forward-backward partons experiments. 
The predicted results for 90°-trigger geometry are similar to those at RHIC. The 



23 



non-trigger Cherenkov gluons should be emitted within the rings at polar angles of 
tens degrees in c.m.s. at LHC by the forward moving partons (and symmetrically 
by the backward ones). This idea is supported by events observed in cosmic rays 
[231 El 121] . The experiments with triggers positioned at various angles to the collision 
axis should be done at the LHC. 

9 Conclusions 

The equations of in-medium gluodynamics (fT5l) . ffl6l) are proposed. They remind the 
in-medium Maxwell equations with non-Abelian terms added. Their lowest order 
classical solutions are similar (up to the trivial color factors) to those of electrody- 
namics (|2"2"|) . especially, for Cherenkov gluons. 

Some effects due to Cherenkov gluons at SPS, RHIC, cosmic rays and LHC 
energies have been discussed. The comparison with experimental data of RHIC 
allows to determine the real and imaginary parts of the nuclear permittivity. It 
is related to the forward scattering hadronic amplitudes. This helps interpret the 
experimental results about the asymmetry of shapes of resonances traversing the 
nuclear medium. Some indications on the energy dependence (dispersion) of the 
nuclear permittivity is obtained from comparison of RHIC and cosmic rays data. 
This asks for the distinction between the different coordinate systems in which the 
Cherenkov radiation (and nuclear permittivity) should be treated for partons moving 
in different directions with different energies. 

Some estimates of properties of the nuclear matter formed in ultrarelativistic 
heavy-ion collisions have been done. This consideration predicts new features at the 
LHC [25]. 
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Figure captions. 



Fig. 1. Away-side azimuthal correlations for STAR in the rapidity and transverse 
momentum intervals indicated in the text and their description by the hypothesis 
about Cherenkov gluons. 

Fig. 2. Away-side azimuthal correlations for PHENIX in the rapidity and trans- 
verse momentum intervals indicated in the text and their description by the hypoth- 
esis about Cherenkov gluons. 

Fig. 3. The distribution of the number of produced particles at different dis- 
tances from the event axis r in the stratospheric event at 10 16 eV [23] has two 
pronounced peaks. Correspondigly, the pseudorapidity distribution possesses two 
such peaks. 

Fig. 4. Excess dilepton mass spectrum in semi-central In(158 AGeV)-In of NA60 
(dots) compared to the in-medium p-meson peak with additional Cherenkov effect 
(dashed line). 
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